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Abstract Stimulated by the recent debate on the 
physical relevance and on the predictivity of g-Gaussian 
formalism, we present specific analytical expressions for 
the parameters characterizing non-Gaussian distribu- 
tions, such as the nonextensive parameter q, expressions 
that we have proposed for different physical systems, an 
important example being plasmas in the stellar cores. 



Keywords Exact results. Rigorous results in statisti- 
cal mechanics. 



1 Introduction 

Many recent experimental measurements of space and 
momentum distribution functions of particle systems 
and of physical quantities in several fields of natural and 
social science show behaviors that can be described with 
non-Gaussian distributions: from cosmic rays to dark 
matter, from astrophysical systems to turbulence and 
quark-gluon plasmas, from biology to finance. More 
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concisely: power laws are detected in several complex 
systems (Gell-Mann and Tsallis 2004). 

Among various non-Gaussian distributions that 
emerge from consistent thcrmodynamical and statis- 
tical frameworks (Kaniadakis et al. 2004; Kaniadakis 
et al. 2005; Lavagno et al. 2007), g-Gaussians, based 
on the so-called nonextensive statistical mechanics in- 
troduced by Tsallis (1988), are appealing for their sim- 
plicity and have found the most applications. A critical 
debate on its theoretical foundation and on its practical 
applicability is still going on. In particular, there has 
been a great effort in this field to construct a statisti- 
cal mechanics capable of describing systems affected by 
nonlocal and memory effects. Relevance of these prob- 
lems in astrophysics and space science has been widely 
discussed in many papers Saxena et al. (2004,2006) 
(see also Mathai and Haubold 2007; Haubold and Ku- 
mar 2007 for comments on applications of g-Gaussians 
in the above field). 

Two important questions regarding q-Gaussian dis- 
tributions have been recently raised by Douxois (2007): 
if g-Gaussian law describes physical phenomena and if 
the entropic q parameter could be predicted in terms of 
microscopic parameters. Tsallis (2008), has answered 
these questions reviewing the present state of nonex- 
tensive statistical mechanics and reporting analytical 
relations for the parameter q that have appeared over 
the recent years. It is our opinion that natural field of 
application and validation of g-Gaussian distributions 
is astrophysics and space science. 

In this work, we present additional interpretations 
of the non-Gaussian parameters, mostly in terms of the 
parameter q, with the corresponding analytical expres- 
sions, that our group have proposed: a general sta- 
tistical interpretation based on Fokker-Plank kinetic 
equations and on Langevin equation, phenomenologi- 
cal interpretations about weakly nonideal plasmas and 
charged particles in electric and magnetic fields inelas- 
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tically interacting with a medium and an interpretation 
based on the quasi-particle description of correlated 
particles with finite lifetime. This list is not exhaus- 
tive, but we believe it offers an important contribution 
to the scientific debate on the physical interpretation 
of the deformation and on the possibility of directly 
linking parameter q, or other analogous parameters, to 
the underlying microphysics; in these specific cases, the 
link has simple analytical forms. In the cases presented, 
the deviation of the distributions from the exponential 
case, q 1, is due to the non-trivial spatial correlations 
among the particles of the system or time correlations 
of long-living metastable systems. The possibility of 
a two-parameter distribution that capture the power- 
law behavior in two different asymptotic regimes (Ka- 
niadakis et al. 2004. 2005) and its application to chaotic 
systems is also recalled (Tonelli ct al. 2006). 

The physically important case of plasmas in stel- 
lar cores is discussed to clarify the importance of non- 
Gaussian distributions. 



2 Non-Gaussian distribution parameters 

By solving a nonlinear Fokker-Planck equation with dy- 
namical and friction coefficients given by a polynomial 
expression with positive coefficients, in 1992 we derived 
(Kaniadakis and Quarati 1993) a set of stationary non- 
Maxwellian distributions with depleted and enhanced 
tails compared to Maxwellian tail. 

Analytical forms of the distributions were given for 
expansion up to the sixth order, with the parameters of 
the distributions expressed in terms of the derivatives 
of the dynamical friction and diffusion coefficients. All 
these distributions deviates from the Maxwellian one, 
which is recovered in the limit of vanishing nonlinear 
corrections to the Fokker-Planck equation. 

As a specific example, Eq. (30b) of Kaniadakis and 
Quarati (1993) 



^ri^i - 1/2) 
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(2.1) 



with ^ = 1 + /3o/(27i), where /3o and 71 are the coef- 
ficients of the leading and subleading term of the ex- 
pansion in the velocity of the dynamical friction and 
diffusion terms. 

Later on, we realized that this specific distribution, 
which in our paper we called second-order stationary so- 
lution of the Fokker-Plank equation, basically coincides 
with the g-distribution obtained independently by Tsal- 
lis a few years earlier with an entropic approach (Tsallis 
1988). In this case the analytical form of corresponding 



entropic parameter q is 

fi 271 + f3o 



(2.2) 



Clearly, non-Gaussian distributions and extended 
entropies were already known and used in different 
fields (Mathai and Rathie 1975; Haubold et al. 2007; 
Kaniadakis and Lissia 2004), however it is only af- 
ter Tsallis seminal paper that systematic extensions of 
thermodynamics and statistical mechanics were devel- 
oped. Such approaches received great attention and are 
actually used in many fields in science. 

In Kaniadakis and Quarati (1997), we have inves- 
tigated more in detail the relationship between our re- 
sults and the nonextensive statistical mechanics results, 
realizing the connection between the physical meaning 
of some of our results of Kaniadakis and Quarati (1993) 
and the corresponding Tsallis distribution and entropic 
parameter q. 

We have also studied Kaniadakis et al. 2004, 2005) a 
wider class of generalized two-parameter statistical me- 
chanics that includes not only the q-thermostatistics, 
but also other physically motivated extensions such as 
the K distribution with its appealing relativistic struc- 
ture (Kaniadakis 2002, 2005), and the Abe distribution 
related to quantum group (Abe 1997). The interpreta- 
tion of the two parameters exist only for selected cases 
or when there remains only one effective parameter and 
it is still object of study. 

2.1 Polynomial expansion of diffusion and drift 
coefficients 

One of the distributions emerging as equilibrium so- 
lutions of the nonlinear Fokker-Plank equation for a 
specific form of the coefficients was first interpreted in 
term of g-distribution in 1997 (Kaniadakis and Quarati 
1997), where we gave an analytical expression of q in 
terms of the ratio of the noise induced drift over total 
drift for specific model with leading order terms in a 
velocity expansion: 



40- 1 



q = 



26- 1 



where 



Jn{v) 



Jniv) + Jdiv) 



(2.3) 



(2.4) 



with Jn{v) = ^ 7 (1 — g) V, J^iv) = ^ {3 — q)v and 7 is 
a constant. The drift coefficient J{v) is the sum of the 
deterministic drift Jd{v) plus the induced-noise drift 

Mv). 

When the noise-induced drift is absent we have q = 
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1 and the q noncxtcnsivc distribution reduces to 
Maxwellian one. The microscopic Langevin equation 
related to Tsalhs distribution is 



dv 



m p L 



and the solution can easily be written in 



1/2 



V(t) 



J r{t)dt + arcsinh (- 



(2.5) 
(2.6) 



where T{t) is a random Gaussian variable with zero 
mean value and S correlation function, v = 2[1 ~ (q — 
1) P X]/[f3m{q-l)],t = 2/[j {q-l)],q = {a+A)/{a+2) 
and A is an energy dimensional normalization constant. 

2.2 Phenomenological derivation for correlated 
particles in plasmas 

Soon after the interpretation in terms of drift and dif- 
fusion coefficients, it was possible to derive an analyt- 
ical expression for q concerning stellar/solar core plas- 
mas, where particles arc correlated and entangled to 
the electromagnetic fields, and their interaction cannot 
be described by local two-body potentials (Coraddu et 
al. 1999). This first specific analytical expression of q 
in terms of microscopic parameters, derived from the 
solution of a Boltzmann equation valid when deforma- 
tion of the global thermodynamic equilibrium distribu- 
tion (Maxwellian) is small (about 0.8 < q < 1.2), con- 
cerns the ionic components of a weakly nonideal neutral 
plasma, as the solar core, in presence of an electric mi- 
crofield distribution. 

The parameter q can be given in terms of plasma 
parameter T ~ (Ze)^n^/'^/(kT), the ratio of the mean 
(Coulomb) potential energy and the mean kinetic (ther- 
mal) energy, where n is the average density, and of a 
parameter a of the ion-ion correlation function (whose 
value is 0.4 < a < 0.9), as defined by Ichimaru (Yan 
and Ichimaru 1986; Ichimaru 1992). By fixing g, we 
can describe with a simple expression the distribution 
of a system of correlated particles interacting elastically 
not only through a pure Coulomb potential but also, for 
instance, through a reinforced Coulomb interaction. 

2.2.1 The effects of random microfields in plasmas. 

The deviation from Gaussian distribution in plasmas 
(Coraddu et al. 1999) has been recently generalized 
(Ferro and Quarati 2005) to consider the presence of 



a microfield distribution of a random force F. The ex- 
pression of q valid for small deformations, 0.8 < q < 1.2 
is 



It 



k /i^ 



(2.7) 



where the two signs indicates subdiffusion or superdiffu- 
sion, k is the energy transfer coefficient, n is the number 
density, fx is the reduced mass, ao and ai are dimen- 
sional constants related to the particle interactions and 
correlations through the cross sections (Jo{v) = ao/v 
and (Ji{v) = ai and therefore to the collision frequen- 
cies vq and vilv). If F is proportional to the electric 
field E {F = e E) and E is greater than a given critical 
value i?crit7 then we obtain 



q = l±2S, valid for (5 < 1/2 , 
where 

5 = 12r^a^ , 



(2.8) 



(2.9) 



with a the above mentioned correlation parameter. 

More concisely, we can derive that the correction S is 
the square of the ratio between the thermal energy den- 
sity over the reinforced energy density times the energy 
transfer coefficient. 

For stellar cores, considering the different ionic com- 
ponents one finds that very small corrections to the 
value g = 1 of the order, let us say, of a few percents 
can sensibly change the thermonuclear fusion rates with 
important consequences on the processes occurring in 
the core and leaving unchanged the bulk properties of 
the star (Coraddu et al. 1999). In the case of the sun, 
helioseismology can test the proton distribution func- 
tion (Degl'Innocenti et al. 1998). 

2.2.2 Electromagnetic fields in two-energy-level 
medium. 

Rossani and Scarfone (2000) have self-consistently de- 
rived from a linear Boltzmann equation an analytical 
expression of the q parameter as an explicit function 
of the electric and magnetic field when charged par- 
ticles are inelastically interacting with a medium en- 
dowed with two energy levels. Concisely 



q^l + -kTf{e, b) 



with 



(2.10) 



(2.11) 



4 



where e — qE/m and b = qB/m arc the rescaled 
electric and magnetic fields, cq and c* are constant co- 
efficients related to the microscopic cross section (due 
to the inelastic contributes), and AE is the energy gap 
between the excited and the fundamental level of the 
medium. 

2.2.3 Particle correlation function and phase space 
cell. 

The expression of the square fractional deviation of 
the number density from a Maxwell-Boltzmann dis- 
tribution of a system of correlated particles by means 
of the evaluation of the phase space volume occupied 
by a nonextensive system of N classical particles that 
slightly deviates from Maxwcll-Boltzmann phase space 
volume gives the possibility to derive an analytical ex- 
pression of q in terms of the radial correlation function 
g{r) (Quarati and Quarati 2003). 

By considering the square fractional deviations from 
Maxwell-Boltzmann distributions defined by V - ^^"'^ , 
where Sui is the variation of the number distribu- 
tion and is the number distribution proportional 
to Boltzmann factor and using the relation (Goodstein 
1975) 



E 



1 



y J [g{r)~l]dr, (2.12) 



where the sum is extended over all energy levels, ne- 
glecting large deviations, we have derived the expres- 



show (Quarati and Scarfone 2007) that, following the 
development usually adopted in superstatistics for the 
inverse temperature(Beck 2001; Wilk and Wlodarczyk 
2007), a neutral astrophysical plasma, as the one of a 
stellar core, has a parameter q defined by 



1 



A(l/i? 



DH) 



1/R 



DH 



(2.16) 



where A(l/i?DH) is the fluctuation of the inverse 
Debye-Hiickcl radius. If we know the equation of state 
of the system the correction factor can be expressed in 
terms of temperature fluctuation. 

2.3 Quasi-particle life-time and non-Gaussian 
momentum distribution 

Many properties of a system of interacting particles can 
often be described by weakly interacting excitations or 
quasi-particles. 
Let us define 



(2.17) 



the self energy of the one-particle propagator, where tu 
and p are related by w = p^/2 m -I- T,{uj, p^). 
When the imaginary part of the self-energy cannot be 
disregarded (Sj > 0) the momentum distribution devi- 
ates from a Maxwellian distribution, even if the energy 
itself is Gaussian. By limiting ourselves to the case of 
g > 1 for small Ej > we have derived (Coraddu et al. 
2000) a phenomenological interpretation of the param- 
eter q as 



(2.13) 



Therefore, when g{r) — > 1 and N ^ oo we have that 
q ^ 1. Otherwise if ^ oo but the particles are 
correlated because we take g{r) ^ 1 we obtain 



9=1-^/ {9(r) - 1] dr , 



(2.14) 



and if the particles arc not correlated \g[r) = 1] but the 
number is small and finite we have 



(2.15) 



2.2. Fluctuations of collective effective parameter. 



q^l + C 



(2.18) 



where C is a constant, Cp — /2 m* + Er and m* = 
m[l- (aER/acj)]/[l + 2m(aER/5p2)]. 

Spatial and temporal correlations among ions have 
large effects on thermonuclear reactions that occur be- 
tween high-energy ions tunneling the Coulomb barrier. 
Quark-gluon plasma is also formed in heavy-ion col- 
lision and similar plasma effects are important in the 
evaluation of reaction rates. This non-Gaussian mo- 
mentum distribution arising from the finite lifetime of 
the quasi-particles should always be taken into account 
when studying nuclear processes in any plasma, the 
more so when the tail of the momentum distribution 
is important. 



Correlations among particles and fluctuations of inten- 
sive parameters, such as the inverse Debye-Hiickel ra- 
dius, l/i?DH, produce nonlinear effects inducing the 
system in a stationary metastablc state that can be de- 
scribed by a g-Gaussian distribution. It is possible to 



3 Two-parameter statistical mechanics: a and 
/3. 

Entropies or power-law distributions of many anoma- 
lous system entropies might be more completely de- 
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scribed by a two-parameter class of logarithms (Ka- 
niadakis ct al. 2004, 2005; ToncUi et al. 2006). In 
fact, the two parameters basically characterize the pos- 
sibly different power-law behaviors of the distribution 
for large positive and negative arguments. 

As an example, at the onset of chaos the logistic map 
shows an universal power-law behavior that fixes one of 
the exponents. This same exponent characterizes the 
time dependence of the entropy and the sensitivity to 
initial conditions: one measures one of the two and can 
predict the behavior of the other. The second exponent 
can be related to the ratio of growth of entropies in 
different formulations (Tonelli et al. 2006). 

This generalized entropies and distributions can be 
introduced by means of a corresponding two-parameter 
class of logarithms 



log(e) 



a + 13 



(3.1) 



and an infinity of one-parameter deformed logarithms 
can be obtained as particular cases by imposing rela- 
tionships between the two parameters. 

In some cases the emerging deformed parameter can 
be related to the microscopic quantities of the system 
under inspection. The g-statistical mechanics is recov- 
ered, as one-parameter subclass of the two-parameter 
statistical mechanics, for a = 1 — g and /? = whilst, 
for instance, for a = /? = k, we obtain a new one- 
parameter subclass of deformed logarithms (Kaniadakis 
2005). Therefore, the question of the analytical deriva- 
tion of the deformed parameter arise also in other con- 
text than the g-statistics. To cite an example, the de- 
formed parameter k has been determined implicitly in 
Kaniadakis (2002. 2005) by considering the black-hole 
entropy, obtaining the relation 



(3.2) 



where Al is the mass of the black-hole. A similar rela- 
tion has also been derived in the case of the (/-statistics. 



4 Application to stellar plasma and 
equilibrium time 

Because many nuclear reactions in the stellar burning 
core proceed by way of quantum penetration of a high 
Coulomb barrier, their cross sections grow exponen- 
tially with energy. Therefore, thermal averages do not 
probe the average energy of the distribution (/cT), but 
its high-energy tail. Consequently, rates are very sensi- 
tive to a relatively small part of the distribution. 



The reasons of such deviations from the exponen- 
tial tail have been discussed elsewhere (Coraddu et al. 
1999; Ferro and Quarati 2005) and include the pres- 
ence of distributions of random microfields larger than 
a critical value and the quantum uncertainty between 
energy and momentum of effective particles which leads 
to a non-exponential momentum distribution even in 
presence of an exponential energy distribution. Ana- 
lytical expressions for these effects have already been 
presented in this paper. 

In addition, screening affects fusion rate in stars not 
only by lowering the Coulomb barrier and, therefore, 
increasing the penetration, but also by modifying the 
distribution. In fact, screening induces effective non- 
local and retarded interactions between colliding par- 
ticles, which can not be treated as pure Coulomb col- 
lisions. In this respect, since the number of particles 
inside the Debyc-Hiickel sphere is too small one should 
also introduce modifications of the usual Debye-Hiickel 
screening. 

Finally, there could be further non-Gaussian effects 
due to the fact that distributions can be out of equi- 
librium. Indeed, nuclear burning in stellar cores is a 
non-equilibrium process with the corresponding ion dis- 
tributions that, in principle, depart from exponential. 
However, the deviation from the pure exponential due 
to the reacting ions is completely negligible for stellar 
structure. As Bahcall et al. (2002) have argued, this 
effect is proportional to the ratio of the collision and 
nuclear reaction times 



T'Coul 
^nucl 



(4.1) 



which is usually very small; for instance, this ratio can 
be estimated of the order of 5 ~ io~ (20-^28) pj-Q^ons 
in the sun. Minor and faster reactions can yield larger 
ratios. 

This line of arguing is based on the assumption that 
only two timescales are relevant for the distribution. 
In fact, collective modes on several scales are present 
in nuclear plasmas that have different timescales. For 
instance, the ratio vi/vco\i\, where ui is the collision 
frequency from the so-called reinforced Coulomb scat- 
tering, is not at all negligible, possibly leading to non- 
Gaussian behavior. 

Therefore, one should carefully assess what is the 
equilibrium relevant to the specific problem at hand 
(nuclear fusion, radiative recombination, diffusion, etc.) 
and also consider, as also argued in Lavagno and 
Quarati (2000, 2006), that non-Gaussian behavior in 
plasmas arises from all collective behaviors discussed 
above not only from non-equilibrium. 

We remind that astrophysical and heavy-ion plas- 
mas require relativistic extensions of the (/-formalism 
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(Lavagno 2002). The reader might be also interested to 
some recent apphcations of noncxtensivc statistical me- 
chanics, where a physical meaning of q is introduced: to 
astrophysics by Jiulin and Yeli (2008a, 2008b) to heavy 
ions by Wilk and collaborators (Osada and Wilk 2007). 

5 Conclusion 

Non-Gaussian distribution functions appear in in- 
teracting particle systems with nonlinear, nonlocal, 
metastable, or long-memory systems; plasmas are typ- 
ical and important example of physical context where 
such distributions appear. 

Modified statistical mechanics, such as the one in- 
troduced by Tsallis, give consistent frameworks that 
naturally yield distributions that deviate from the 
Maxwellian. For many purposes a single effective pa- 
rameter, such as the entropic parameter g, is sufficient 
to characterize the non-Gaussian behavior of the sys- 
tem. This parameter can always, in principle, be de- 
rived from the underlying microscopical dynamics, even 
if complete dynamic calculations are very difficult for 
most of the systems. 

Even when the calculation from first principle is pos- 
sible, the answer does not need to be a simple analytical 
formula, but one often can find numerically the param- 
eter(s) that describe(s) the distribution; nevertheless 
the result is predictive, since the same number applies 
to many distributions for similar systems. 

In this paper, we have reported a few cases where 
there exist an analytical formula for the entropic pa- 
rameter q in terms of properties of the systems. All 
these cases involve particles that are properly corre- 
lated. Other approaches are certainly possible, but 
nonextensive q statistical mechanics appears particu- 
larly well-suited to describe such complex systems, par- 
ticularly astrophysical and space systems. 

At last, we recall that g-Gaussian distributions, in 
addition to provide a very powerful effective descrip- 
tion of complex systems, have also a predictive power. 
For instance, Lutz (2003) proposed an experiment on 
optical lattices; such an experiment was later performed 
(Douglas et al. 2006) confirming the expectation. We 
have also proposed (Macro et al. 2006) an experiment 
to look for a cut-off in the energy distribution of gam- 
mas in electron recombination that could confirm de- 
viations and possibly measure the relevant parameter 
of the distribution. To obtain predictive results in the 
astrophysics and space science field is the challenge of 
the next future. 
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